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The topological approach to baryon-antibaryon production in the chiral phase transition is nu- 
merically simulated for rapidly expanding hadronic systems. For that purpose the dynamics of the 
effective chiral field is implemented on a space - rapidity lattice. The essential features of evolu- 
tions from initial 'hot' configurations into final ensembles of (anti-)baryons embedded in the chiral 
condensate are studied in proper time of comoving frames. Typical times for onset and completion 
of the roll-down and exponents for the growth of correlations are extracted. Meson and baryon- 
antibaryon yields are estimated. For standard assumptions about initial coherence lengths they are 
compatible with experimental results. 
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dc I. INTRODUCTION 

<: 

7-H ' The topological approach to baryon structure and dynamics in the framework of an effective action for mesonic 
£NJ , chiral fields has achieved a number of remarkable successes. The soliton concept [lj for baryons provides an impressive 
account of spectrum and properties of baryon resonances (essentially without numerous 'missing resonances') [2J, with 
predictive power that recently has even led to the first indications for pentaquarks Model- independent relations 
between T-matrix elements for meson-baryon scattering Q and explicit results for specific channels are well supported 

j. j by experimental data [f| . The matrix element of the axial singlet current related to the spin content of the proton is 

I . naturally of the observed order of magnitude |6(. The 'unexpected' behaviour recently found |7j in the ratio of electric 
and magnetic proton form factors has been predicted in this approach long ago @. The underlying chiral effective 
action is profoundly based on the l/iV c -expansion of QCD preserving all relevant symmetries. Efforts to include 
next-to- leading order quantum corrections have brought substantial improvement as expected |lOj . 

The manifestations of a chiral phase transition pose another natural challenge for an effectice theory with a ground 
state that is characterized by spontaneously broken symmetry. The possible formation of disoriented domains [TlT ] 
CIh! during the growth of the chiral condensate has been in the focus of interest for some time. But signatures in terms of 
anomalous multiplicity ratios for differently charged pions have not been observed [l3 | , in accordance with theoretical 
conclusions |!3lll4|. Anomalies in anti-baryon production were very early recognized as possible signals for interesting 
dynamics |l5j | in that phase transition, and the concept to consider baryons as topological solitons in a chiral condensate 
should lead to quite definite expectations for this process. 

Meanwhile, in relativistic heavy-ion collisions at RHIC, very high energy densities are being produced in extended 
spatial regions which are essentially baryon free and well separated in rapidity from the nuclear slabs receding from 
the collision volume. The experimental values found in the central rapidity region for the ratio of the integrated p to 
7r~ yields lies between 0.065 and 0.075 0. This is still too close to the thermal equilibrium P/tt~ ratio (for a typical 
plasma temperature of T ~ 200 MeV) , 

P/tt- ~ 2exp((m ;r - m p )/T) = 0.035 (1) 

to constitute a clear indication for interesting underlying physics. Still, although the experimental result does not 
look very exciting, it still poses a constraint for the possible validity of the soliton concept, because any conceivable 
dynamical production process must be able to produce a comparable number. 

In the topological approach the number of baryon-antibaryon pairs produced during the chiral phase transition 
depends on two factors: the first is the modulus \p\ of the average winding density present in the initial 'hot' field 
configuration. In analogy to applications in cosmology fr7j| and condensed matter systems |l8j this quantity is closely 
related to the coherence length for the local orientations of the chiral field 3>. Without detailed knowledge about the 
initial field configurations this coherence length enters as a parameter and takes away stringent predictive power from 
the approach. However, different conjectures about the nature of the initial field ensemble suggest typical ranges for 
the coherence lengths which then may be discriminated by the experimentally observed abundancies. 
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The second factor is the reduction of the initially present total n>i = J \p\dV through the dynamical ordering process, 
which finally leads to the formation of stable soliton structures embedded in the topologically trivial ordered chiral 
condensate of the 'cold' system. The solitons or antisolitons evolve from topological obstacles which are met by the 
aligning field orientations, and develop into their stable 'cold' form during the course of the evolution. At the end, the 
same integral rif — J \p\dV counts the number of finally surviving nontrivial separate structures, so it is identified with 
the number of baryons and antibaryons created in the process. The decrease of n during the roll-down is reasonably 
well represented by a power law (t/to)~ 7 and the exponent 7 can be measured in numerical simulations. Evidently, 
the initial time tq which marks the onset of the evolution, enters here as a second parameter which further reduces 
the predictive power of the approach. Fortunately, it turns out that 7 is rather small, so the dependence on tq is only 
weak. 

Measuring 7 and the time Tf when the roll-down is completed, presents a typical task for numerical simulations once 
the equation of motion (EOM) which governs the field evolutions is implemented on a lattice. The underlying effectice 
chiral action is known from other applications, so no additional parameters enter at this point. In condensed matter 
applications, a phase transition is generally driven by an externally imposed quench, or by a dissipative term included 
in the EOM. In cosmology or in our present heavy-ion application it is the rapid expansion of the hot volume which 
drives the cooling process. This expansion is efficiently implemented [l9| by transforming to rapidity - proper-time 
coordinates, i.e. by boosting to the local comoving frame. This is especially convenient if we consider a system that 
expands only in one (longitudinal) direction with its transverse scales unchanged, (the Bjorken rod). The resulting 
dilution of the longitudinal gradients drives the system towards its global minimum. However, as there is no genuine 
dissipation in the system, the total energy approaches a constant which resides in the chiral fluctuations around the 
global minimum. Thus, the simulations also allow to estimate pion- or sigma- meson abundancies. 

Naturally, before the field configurations can roll down towards the global minimum, the potential V(§ 2 ,T) which 
underlies the EOM must have changed from the 'hot' chirally symmetric form to its 'cold' symmetry- violating form. 
But, during the early stages, the evolutions are dominated by local aligning of the field orientation <&. During this 
phase the form of the potential is not important. So its time dependence can be replaced by a sudden quench where 
the 'hot' field configuration is exposed to the 'cold' potential V(§ 2 ,T — 0), from the outset at initial time To . In the 
following, for definitcness we make use of this sudden quench approximation, (although the simulations, of course, 
allow to study other cases as well). 

For the sake of simplicity we first discuss all relevant features for the case of the 2-dimensional 0(3)-model, with only 
one spatial dimension transverse to the longitudinal rapidity coordinate. Except for computational complexity the 
extension to the 3-dimensional (3(4)-field presents no essential new features. The effective action, its transformation 
to the Bjorken frame, and the resulting EOM are presented in section II. It is important for the choice of the initial 
ensemble of field configurations that it allows in a convenient way to monitor the initial coherence lengths because they 
are the crucial parameters for the final baryon-antibaryon multiplicities. We choose an isotropic Gaussian random 
ensemble of field fluctuations in momentum space which is characterized by a temperature-like parameter to be able 
to compare with other approaches. Of course, this is not necessary. In fact, even at initial time tq the longitudinally 
expanding Bjorken rod need not be an isotropic system, and it may be physically justified to distinguish already in 
the initial ensemble two different, longitudinal and transverse correlation lengths. This is easy to incorporate, but in 
section III we present initial conditions which are locally isotropic. 

As discussed elsewhere pfl | stable solitons shrink in a spatially expanding frame. Therefore, lattice implementations 
of their dynamics will necessarily involve lattice artifacts after some time. These are discussed in section IV. They 
can be isolated and subtracted from the physically interesting quantities. 

In section V the essential features of typical evolutions are discussed. Estimates for the times of onset and completion 
of the roll-down are obtained, and the dynamical exponents for the growth of correlation lengths and decrease of defect 
number are established and compared. The spectrum of the fluctuations remaining after the roll-down is considered 
and finally the mesonic and baryonic multiplicities are obtained. 

The extension to the physically interesting 3+1-dimensional (3(4)-field is discussed in section VI. The topological 
generalization is well known, the additional transverse dimension is of little influence for the growth exponents. 
However, the coupling constants in the effective action here are related to physical quantities, so they are known 
(except for some uncertainty concerning the cr-mass) , and the results can be compared with experimentally determined 
abundance ratios. 

Of course, it would be desirable to obtain a very definite answer whether the topological approach to antibaryon 
production in a chiral phase transition is validated or ruled out by the data. However, with our poor knowledge about 
the initial conditions in the hot plasma after a heavy-ion collision, we cannot expect much more than allowed ranges 
for the relevant parameters, which hopefully overlap with standard ideas about coherence lengths and formation times. 
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II. THE EFFECTIVE ACTION IN THE BJORKEN FRAME. 

For simplicity we first discuss the 2+1 dimensional 0(3) model. It is defined in terms of the dimensionless 3- 
component field <I> = <3><1> with unit- vector field 4>, • $ = 1), and modulus ('bag'-)field $, with the following 
lagrangian density in 2 + 1 dimensions [x, z, t) 

C = fi(cW+CW+C<®) (2) 



(/ 2 is an overall constant of dimension [mass 1 ], so the physical fields are of mass-dimension [mass 1 / 2 ]). The 
second-order part comprises the kinetic terms of the linear a model 

£< 2) = ^Stf**, (3) 

£( 4 ) is the four-derivative 'Skyrme'-term (which involves only the unit- vector field €>) defined in terms of the topological 
current p^ 



(? = -U^* ■ (d v & x d p $), (4) 

07T 



(which satisfies d^p^ = 0) 

£( 4 ) = _A£ 2 Qli( f 



2A/ 2 r „ ~ 

"" {d^d^&f - (d^d^^^d"®) , (5) 



(87T) 2 

and C^°> contains the <I> 4 potential and an explicit symmetry-breaker in 3-direction 



= -V(*,T) = -~ (* 2 - /(T) 2 ) 2 - if$ 3 ) - const. (6) 



^ 2 I 4 

with dimensionless coupling constants A and H, and 

/ 2 (T) = / 2 (T)--^— . (7) 

This choice ensures that the global minimum of the potential V(&,T) is always located at <&o = (0,0, /q(T)). 
Generically, the function /g(T) decreases from /q = 1 at T = towards zero for large T. The constant in the 
potential I© is chosen such that the value of the potential V at $ = is independent of T, (given by the constant 
V(0,T) = (A + 2iJ)/(4^ 2 )), and at the (T = 0)-minimum * = * = (0,0,1) we have V{$ ,T = 0) = 0. 
The masses of the w- and cr-fluctuations (tti , Hi , fa + a) around this minimum are 

2 H 2 ^A/q 2 /„<. 

Without explicit symmetry breaking, H = 0, we assume that f 2 (T) changes sign at T = T c , such that <&o = (0,0,0) 
and m 2 = m 2 = to 2 = A|/ 2 |/^ 2 for T > T c . 

The parameter ^ (with dimension of a length) which we have separated out from the coupling constants of potential 
and Skyrme term can be absorbed into the spatial coordinates x. So it characterizes the spatial radius of stable 
extended solutions (which scales like l/\f~p)- As i simply sets the spatial scale, it could be put equal to one, as 
long as no other (physical or artificial) length scales are relevant. For lattice implementations, however, the lattice 
constant a and the size of the lattice (Na) set (usually unphysical) scales. To avoid artificial scaling violations we 
have to ensure that the size of physical structures (like solitons) is large as compared to the lattice constant a and 
small as compared to the lattice size Na. So, for numerical simulations we have to choose 1 <C t/a -C N. It has been 
shown in ref.[21j that for solitons which extend over more than at least 4-5 lattice units the energy Eb is independent 
of £/a. So, in the following we will adopt £/a ~ 5 as sufficiently large. This appears also as physically reasonable, if 
we consider typical lattice constants of 0.2 fm and baryon radii of about 1 fin. On the other hand this will require 
lattice sizes of at least N ~ 50 to avoid boundary effects for the structure of individual solitons. Unfortunately, in 
the Bjorken frame which we shall use in the following, the longitudinal extension of stable solitons shrinks like r _1 
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as function of proper time r. This means that after times of order I the simulations will be influenced by lattice 
artifacts, which may even dominate for large times. 

For rapid expansion in (longitudinal) z-direction we perform the transformation from (z, t) to locally comoving 
frames (t],t) with proper time r and rapidity rj, defined through 



t 



rcosh?7, 
rsinh?7. 



<1 



y/t 2 - z 2 , 
atanh y— ^ 



d t = cosh 77 d T — S ^ n ^ ^ Q 

T 

d z = — sinh^9 r H — -d v 



(9) 



Inserting © into and (JSJ leaves the form of and invariant, with <9 t replaced by d T , and d z replaced 
by -d v . The specific structure of the Skyrme term again eliminates all terms with four r- or ^-derivatives. For the 
effective action we take the integration boundaries from — oo to +00 for rapidity 77 and for the transverse coordinate 
x. The 3-dimensional space-time volume element dx dz dt is replaced by r dx drj dr. Therefore, in a separation of the 
action S in kinetic terms T, gradient terms L, and the potential U, 



S 



dT 



+ OO 



£dr/dx= / (T± + Tii -L±- L11 - U)dr 



(10) 



the longitudinal |j-terms involving rapidity gradients carry a factor 1/r, while all other terms carry a factor r. So we 
have 



Tii = - 




(<9 T $ x d x $) 



drj dx, 



-1 2 ~ 



-p • (a T $ x 



in 



(4tt) 2 



A .^2 ^ 



cZ?7 <ix, 




$3 



drj dx, 



4^2 (* -/ 2 ) -72 ^3 + const.}. df7dz 



f 2 



(11) 

(12) 
(13) 
(14) 
(15) 



Variation of 5 with respect to $ leads to the equation of motion (EOM) . 

The contributions of D- 4 ' to the longitudinal and transverse parts Ty and T± of the kinetic energy cause certain 
numerical difficulties for the implementation of the EOM on a lattice. They require at every timestep the inversion 
of matrices which depend on gradients of the unit- vectors 3>, which multiply first and second time-derivatives of the 
chiral field. This can be troublesome in areas where the unit-vectors are aligned, and can be poorly defined in regions 
where the unit- vectors vary almost randomly for next-neighbour lattice points (i.e. for initially random configurations, 
or near the center of defects). In any case, stabilizing the evolutions requires extremely small timesteps and leads 
to very time-consuming procedures. Although these problems can be handled, we have compared the results with 
evolutions where the kinetic energy is taken from alone. For coupling strenghths X£ 2 within reasonable limits, 
we find that the resulting differences do not justify the large additional expense caused by the fourth-order kinetic 
contributions. Evidently, the reason is, that the EOM determines the field-velocities (depending on the functional 
form of the kinetic energy) in such a way that the numerical value of the total kinetic energy is not very sensitive to 
its functional form. We therefore use in the following an effective action where the kinetic terms <|11H and (|12fl are 
replaced by 



(9 r $3 r $) drjdx, 



With this simplification the EOM is 
1 



-<9 T * + <9 rT * - 9 IX * 



1 



2 d nv & 



X 



Tu = 0. 



H 



(16) 



^2 



0. 



(17) 



This form has the big advantage that we can make use of the geometrical meaning of the winding density po as the 
area of a spherical triangle, bounded by three geodesies on a 2-dimensional spherical surface. In closed form it is 
expressed through the unit-vectors pointing to its corners, and does not involve gradients. So this allows for a very 
accurate and fast lattice implementation of the last term in the EOM. 
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III. INITIAL CONFIGURATIONS 



We assume that at an initial proper time To the system consists of a hadronic fireball with energy density £o stored 
in a random ensemble of hadronic field fluctuations. Subsequently, for r > r , it is subject to EOM (|17|l . The 
initial condition and the symmetry of the action imply boost invariance, i.e. the system looks the same in all locally 
comoving frames, so it is sufficient to consider its dynamics in a rapidity slice of size Ar/ near midrapidity r/ = 0, which 
constitutes a section of the initially created Bjorken rod with transverse extension A. The energy E = T + L + U in 
this slice then is given by an ^-integral which extends over the finite rapidity interval A77 and represents the energy 
contained in a comoving volume V = rArjA. Due to the symmetry of the initial condition this comoving volume 
grows with increasing proper time t into spatial regions with high energy density, therefore E contains contributions 
which increase with r. The average energy density e = E/V satisfies de/dr < 0. 

For numerical simulations we implement the configurations &(x,r),T) on a rectangular lattice (x,i]) — (ia,jb) 
— 1...N) with lattice constants a for the transverse coordinate and b for the rapidity lattice. We define the initial 
configurations <&y at the lattice sites {i, j) as Fourier transforms of configurations 3?^ on a momentum lattice 



f 



N/2 



*V = jf E J(e i * W, *H+".), (18) 

k,l=-N/2+l 

with <& fc j = <fr_fc_2. Inversely, the real parts otki and the imaginary parts (3 kl of &ki are obtained from the real 
configuration through 



1 - . 2ir 

a ki = E C0S "/V ^ k + * y = a - k - 1 ( 19 ) 

ij = l 

1 N 2tt 

Pkt = ~Jr E sin — (ik + jl) * y = -(3_ k _ t , (20) 
so we obtain the spectral power P pq of the configurations (or a specific component of it) at any time r from 



P pq = * fc ; • = (Xki ■ ol u + /3 kl ■ /3 kl (21) 

for any transverse or longitudinal momentum (p, q) = ^jr(k, I), for (k, I = —N/2 + 1, N/2). 

For the initial configurations at r = To the real and imaginary parts of each of the three components of <&ki at each 
momentum-lattice point (p, q) are chosen randomly from a Gaussian deviate Gfcz($) with /cZ-dependent width Uki, 

G w (*) = -^=ex P (-^-|, with <7& = ^<>>q>| -L _). (2 2, 




with normalization Z chosen in such a way that 

N/2 

E < = N 2 *l (23) 

k,l=-N/2+l 

(In the continuum limit (a -> 0, N -> 00) we have Z=^(l + f) e" m / r .) 

In other words, we choose a Boltzmann distribution for the average occupation numbers nu = ((^ki^li)) = &ki f° r 
each field component, as for relativistic (non-interacting) particles with mass m. Here the mass m 2 is defined by the 
absolute value 

m 2 (T) = ±\f(T)\ (24) 

for the fluctuations around <& = in the symmetric potential (jSJ) at the initially high temperature T = 7o, where 
/ 2 (T) is negative. The amplitude a 2 , pl a Y s the role of a fugacity 



°o = exp(-M/^") 



(25) 
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for negative chemical potential \x. In the temperature range which we consider (0.05 < aT < 0.8) (cf. fig-®) a 
suitable value for /x is a/x ~ —0.6. (With this choice the average amplitude of the chiral field is not subject to abrupt 
deviations from its initial value immediately after the onset of the dynamical evolution). 

We assume isotropy of the initial ensemble with respect to rotations in 0(3)-space such that the three components 



of the field fluctuations (a = 1,2,3) have the same average square amplitude a^. By picking each component 



independently at each point (k, I) from the Gaussian ensemble, different components are uncorrelated and equal 
components at different points (on the momentum lattice) are also uncorrelated, 



S-kk'S- 



a kl^al3 ( -(Skk'StV 



«') + -^(Skk'Sw - 5- 



kk' 



iS-w) = cr k i5 a fjSkk'Su 



Together with (|18|l this leads to the fluctuation in the real field configurations 



(26) 



i 



N/2 
kd=-N/2+l 



= 6, 



«/3 ^0 



(27) 



which is, of course, independent of the lattice point Its magnitude is controlled by the constant cr 2 , m (EU- It 

should be sufficiently small to keep the amplitudes of the average initial fluctuations small. On the lattice the upper 
limit for the momenta p,q is — , (i.e. fc, I = N/2). So, as long as 



r « 



(28) 



the lattice cut-off (upper limit momentum) imposed by the finite lattice constant is unimportant because the cor- 
responding states are almost unoccupied. Note that periodicity and antisymmetry of the imaginary parts in 1I9|) 
requires that (3 kl vanishes if both k and I are multiples of N/2. With the condition (|28|l satisfied, this holds with good 
accuracy also for the initial configuration picked randomly from the ensemble l|22|l . 

The average number of topological defects in a random ensemble of vector configurations is closely related to the 
characteristic angular coherence length in that ensemble. Therefore, it will be necessary to measure the (equal-time) 
correlation functions for the unit- vector fields <fr for the evolving ensembles. In order to have an analytical result 
at least for the initial configurations (where length and orientation of the 3- vectors are uncorrelated), it is easier 
to consider the correlations among the full vectors <I>. Therefore, we define normalized transverse and longitudinal 
correlation functions 



C±(i) = 



C\\(i) 



I 



3(7^2 



3<rgiV 2 



N N N 



m,n— 1 
N 



m,n— 1 
N 



k,l=l 
N 



k,l=l 



(29) 



with transverse coherence lengths R± and longitudinal (dimcnsionless) coherence rapidity Rii defined through 



C(i)< 



1 



for i > 



R± 



or i > 



(30) 



respectively. For the initial ensemble 122(1 the correlations are, of course, isotropic on the lattice, i.e. 



R± _ R\\ _ Ro 
aba 



(31) 



with initial spatial coherence length Rq. In the continuum limit (a — > 0, N — + oo), we obtain C(r) as function of the 
spatial distance r (or rapidity rj = r(b/a)) 



C(r) 



( 1 + r 2 T 2)3/2 



T 

1 + f 



(32) 
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Specifically, putting m = 0, the coherence length R as defined in (|30|) is 



R=^±~°-f. (33, 

This allows to put limits on the range of temperatures which can be reasonably represented on the lattice. Typically, 
for lattice size of N ~ 100, T should lie within the range from about 0.02 to about 0.8 inverse lattice units. For 
smaller values the inital coherence length already covers more than half of the lattice so almost no defects will fit on 
the lattice, for larger values the correlation lengths approach the lattice constant. It may be noted that with (|24(l . for 
(£/a) ~ 5 and (aT) ~ 0.1, the ratio m/T is not very small, so generally we expect appreciable deviations from the 
T- 1 scaling in 0, (e.g. for (i/a) = 4 we find (R/a) ~ (aT)" 8 , cf. figEJ. 

During the evolution in the Bjorken frame the correlations rapidly become anisotropic. We then conveniently define 
an average coherence length R through 




R 2 

This may be compared to the coherence radius obtained from the angular-averaged correlation function 

1 



C(r) 



N N N 

m,n— 1 i,j m,n— 1 k,l — 1 



(34) 



(35) 



where the i, j-sum indicates an average over all lattice points in a narrow circular ring with radius r around the lattice 
point mn. 

The essential characteristics of the evolutions are not very sensitive to the choice of the initial time derivatives. 
(They can as well be put to zero.) The equations of motion very quickly establish appropriate velocities. Of course, 
the absolute value of the total energy depends on that choice. For the simulations presented in the following we 
construct in analogy to the initial configurations l|18|l an initial ensemble of time derivatives through 

(flr*)« = ^ £ ^(e [ ^ k+ ^* k i-c.c), with u kl = J(^k)\(^l) 2 +rn*. (36) 

The Fourier coefficients &ki again are picked randomly from the same Gaussian deviate l|22l) . 



IV. SHRINKING SOLITONS IN COMOVING FRAMES 

Let $^(x, z) be a static soliton solution of the model in its (x, z) rest frame, which minimizes the static energy 
E = L + U with a finite value for the soliton energy E = E n . After the transformation to the Bjorken frame, the 
configuration ^^'(x, rf) = ^ s \x,rr]) then describes a static solution of the action in the comoving (x,r]) frame at 
proper time r (where d z is replaced by (l/r)9„), for the same value of E . It represents a soliton with the same 
finite radius in transverse x-direction as before, but with its radius in longitudinal ?/-direction shrinking like 1/r with 
increasing proper time r. The total energy Eq of this shrinking soliton is, of course, independent of r. (Naturally, 
this consideration strictly applies only to the adiabatic case, where r is considered as a parameter. In the dynamical 
ordering process the evolution of the solitons towards their static form may appreciably lag behind the actual progress 
of proper time.) For lattice implementations, with the typical spatial radius of the stable solitons given by £, the 
longitudinal extension of the solitons for times t » f has shrunk down to (dimcnsionlcss rapidity-)lattice-unit size 
and longitudinally adjacent solitons no longer interact. In transverse direction, however, the solitons develop their 
stable size of Ija lattice units, they keep interacting, attracting close neighbours or annihilating with overlapping 
antisolitons (cf. fig^]). 

For solitons shrinking longitudinally down to lattice-unit size the energy will begin to deviate from the value Eq as 
soon as the longitudinal extent covers merely a few lattice units. To get an approximate idea for the energy limit let 
us assume that a single separate soliton finally degenerates into a transverse string of 21 + 1 lattice points, on which 
|<I>| varies from nearly zero (in its center) to the surrounding vacuum configuration $o = (0, 0, 1), i.e. <I> = (0, 0, \i\/t) 
for — I < i < i , on that string of lattice points. Then we find for the contributions of a single soliton 

tf~Al. (37) 
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FIG. 1: Soliton configuration after a typical evolution on a 50x50 lattice (for A = l,£/a = 4, <ro = 0.2, H = 0.2, aT = 0.2) 
at time t/to = 1000, i.e. long after completion of the roll-down. The bag-field |#| of the solitons (upper part) is sqeezed 
longitudinally to lattice-unit size; the positive or negative winding densities (lower part) are located at the center of the bags. 



So, apparently, solitons shrinking on a lattice contribute to the energy terms which rise linearly with proper time r 
which (as lattice artifact) will dominate the total energy for large r. 

We expect the winding density of the squeezed defect to be located on v lattice squares near its center. This implies 
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for the fourth-order term 

= (38) 

ii T 

The winding density is determined by the orientation of the field unit-vectors alone, so it sufficient to consider the 
unit-vectors <i>. We expect the squeezed defect to consist of just one unit-vector $ = (0, 0, —1) at the soliton center 
looking into the direction opposite to all surrounding unit-vectors (0, 0, 1). That lattice point is the top of four adjacent 
rectangular triangles (with the diagonals connecting the four nearest neighbour points as bases) which together cover 
an area of two lattice squares. So we expect a winding density p=l/v with v ~ 2. 




FIG. 2: Potential energy U, kinetic energy T^ 2 \ transverse and longitudinal (second-order) gradient terms Lj 2 ' and , 
the number of defects n, and the average length of the chiral field |4>|, for a typical evolution after a sudden quench (for 
A = l,l/a = 4, (Jo = 0.2, H — 0.2, aT = 0.2,7V = 50). For comparison, the straight lines given by eas. l|45^ with !4Ht are 
included. 



This dominance of lattice artifacts for r 3> I is illustrated in fig. © which shows a typical evolution on a 50 x 50 
lattice for I = 4. The total winding number is B = —1. After the roll-down the number of solitons stabilizes 

at n = 9 (cf. figCJ. Apparently, both U and approach a linearly rising limit for r ^> 10 2 , approximately like 
~ \n{rji) which dominates the total energy, but does not affect the (essentially constant) kinetic energy. Longitudinal 
contributions drop off like r _1 , so they are irrelevant. 

It appears from fig.© that for this evolution the roll-down (where the average of $ aproaches the vacuum value 
$ = 1) takes place during the time interval 2£ < t < il, i.e. long before artificial lattice effects dominate the energy. 
It is also by the end of the roll-down that the number of created defects stabilizes. So we would conclude that results 
obtained from lattice simulations for baryon-antibaryon production during the chiral phase transition in a rapidly 
expanding chiral gas are not severely affected by lattice artifacts. On the other hand, to follow the evolutions beyond 
the end of the roll-down, which comprise small 'a' and V-oscillations of $ around the true vacuum, interfering with 
small oscillations of the bag profiles (resonances), will require to subtract the lattice artifacts. 
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V. EVOLUTION UNTIL FREEZE-OUT 



In this chapter we will follow typical evolutions of the chiral field after a sudden quench in more detail and try to 
analyze their characteristic features up to the end of the roll-down. 

Immediately before the sudden quench at r = To the initial ensemble is prepared as described in section II. The 
average length of one component of the chiral field is given by Co (cf. eq. J37J|), the potential in JBJ) is characterized 
by a negative value of f 2 . So, for sufficiently small Cq we have at t = r 



U ° = 7 »-&i / (/ 4 + 2|/ 2 |«* 2 )))dx 



d v = (Co + C 2 )V , (39) 



where Vq = tqAt/A is the inital volume of the Bjorken slice, and the constants are 



For the derivatives at lattice points Eq. i(T%|) implies 

N/2 . . . N/2 . . 

= £ E (w) em ^'^ u ' = £ E (41) 

k,l = -N/2+l ^ ' k,l=-N/2+l ^ ' 

Again replacing the integrands in i|13|l . (|14fl by ensemble averages leads to the second-order gradients contribution at 

T = T 

4 2) = ^ (2) V . (42) 
For the constant C^ 1 we have in the continuum limit 

CW = r« T2 ( 1+ 3T^W))- <«> 

Similarly, one may obtain a rough estimate for averaged over the initial ensemble by replacing in the unitvectors 
•I by <J>/cr - 

During the very early phase of an evolution in proper time the initially random ensemble of fluctuations will 
essentially stay random. This means that the integrals in (| 1 3|l — (| L 5|) will remain constant, given by their initial values. 
Therefore, the time dependence of the different contributions (f 1 3|1 - (f 1 to the total energy is given by the kinematical 
factors (t/to) or (to/t) alone, with the integrals approximated by replacing the integrands through their averages in 
the initial ensemble. 

After the quench, f 2 is positive, so for sufficiently small a 2 we have 



(44) 



and 



L f) = L C ( 2 ) V , = H V . (45) 

T " r 



In fig.© both straight lines, l|45|l with 143|) . are included for comparison. It may be observed that the integral 

(2) 

L| ; involving the longitudinal gradients follows the straight line decrease almost until the onset of the roll-down. 

(2) 

This means that the rapidity gradients basically stay random. On the other hand, the integral L y ± ' follows the linear 
rise only for about one unit of proper time after the onset of the evolution. Already near t/tq ~ 2, the transverse 
gradients are strongly affected by the dynamics and interfere with the kinetic energy. Due to the relative factor of 
1/t 2 of Lit as compared to L± the dynamics quickly gets dominated by the transverse gradients alone, such that the 

(2} (2) 

average kinetic energy follows the average transverse-gradient energy L ± , while the rapidity gradients (in and 
L^) which decrease like 1/r are no longer relevant for the overall dynamical evolution. 
Disregarding rapidity gradients altogether, the EOM l|17|) reduces to 

-a T * + a Tr *-a ra *-m 2 * = o, (46) 

T 
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which describes wave propagation in transverse direction, A(t) e~xp(ipx). Here the mass m 2 again characterizes the 
fluctuations around <I> = 0, 

m 2 = Xf 2 /£ 2 (47) 

so m 2 is negative for negative f 2 (where $ = is the stable minimum), and it is positive for potentials which actually 
do have a lower symmetry-breaking minimum. The amplitudes A(t) generically are Bessel functions, 



A(t) ~ Jo{ry/p 2 - m 2 ) for p 2 - m 2 > 0, 

A(t) ~ hij\Jm 2 -p 2 ) for p 2 -m 2 <0. (48) 

For large values of their arguments the amplitudes of Jo decrease like 1/\/t, while Iq contains exponentially rising 

(2) 

parts. Modes with large transverse wave numbers contribute most to L ± ■ Therefore, with their amplitudes decreasing 

(2) (2) 

like 1/y/r, the kinematical factor r in L ± is compensated. So we expect that the linear rise of L ± ends as soon as 
the dynamics is dominated by the transverse gradients and is followed by a phase where 

«T))^«L ( L 2) ))~const.| T . (49) 

For negative m 2 no amplification occurs. After the quench, however, when f 2 has become positive, a few modes with 
small transverse wave numbers will start to get amplified. Typically, for wave numbers p — 2nk/N, with k integer 
(0 < k < N/2), waves with k/N < \f\f /(2irt) get amplified, e.g. the lowest three or four out of N = 100 for I ~ 5 
(for A = 1 and f 2 — 1). At first, the rate of amplification is slow because the exponential rise is compensated by a 

(2) 

decreasing function for small arguments in Iq(x). These low- A: modes do not contribute much to L]_ . In fact, the 
k = mode, which experiences the largest rate of amplification, does not contribute at all. 

(2) 

While the amplification effect is not very pronounced for L ± , the few slowly exponentially rising contributions 
from the lowest-momentum transverse waves cause a noticeable rise of the condensate (($ 2 )) after some time. This 
enters into the fluctuating part C2 of the potential U and drives it away from its linear rise given by (|44|l . Then also 
the fourth-order terms in the potential become important and the dynamical evolution subsequently is dominated by 
the local potential. This initiates the roll-down of the field configuration at the majority of the lattice points into the 
true vacuum 3>o = (0,0, 1). The transition into the symmetry-violating configuration takes place, with formation of 
bags and solitons in those regions where the winding density happens to be high. 

To estimate the time t\ for the onset of the roll-down we consider the k = mode with amplitude I^irm). 
Amplification of this amplitude by a factor e in the time interval from to to t± requires 

lnJofam) = 1 + lnI (T m). (50) 

The r.h.s. depends only very weakly on To, as long as {Tom) < 1. In fact, (rim) varies only from 2.26 to 2.55 for 
< ("Tom) < 1. So, for convenience we simply take (rim) « 5/2 if (rom) is of the order of 1 or less. Otherwise, for 
larger values of (rom), T\ has to be obtained more accurately from l|50l) . Typically, therefore, the transition from the 
gradient-dominated to the potential-dominated phase, happens near 

(51) 




However, up to this time T\ of the onset of the roll-down, i.e. throughout the whole gradient-dominated phase the 
potential plays no significant role. The overall evolution proceeds practically independently from the (positive or 
negative) value of f 2 in the <I> 4 -potential ©. This also implies that the quenchtime (the timescale for changes in 
f 2 ) is irrelevant as long as it is smaller than the time during which the gradient terms dominate the evolution and it 
justifies the use of the sudden quench approximation where we impose the 'cold' (T = 0) potential from the outset at 
r > t . 

With £/tq > 1, the ratio (ti/to) 2 is sufficiently large to render all longitudinal (rapidity) gradients unimportant 
as compared to the potential. This means that during the subsequent roll-down different rapidity slices become 
effectively decoupled, and begin to evolve independently from each other, while in longitudinal directions the solitons 
contract to lattice unit size. Within these rapidity slices, tt- and cr-modes propagate transversely, and eventual further 
annihilations of soliton-antisoliton pairs take place while the transverse shapes of the sqeezed bags are established. 

By the end of the roll-down the remaining nontrivial and sufficiently separate structures have essentially reached 
their stable form. Apart from small fluctuations, the integral (or sum) over the absolute values of the winding density 
n = J\p\ dxdz then stabilizes and counts the (integer) number of these defects. Therefore we identify the end-of-the- 
roll-down time with the (chemical) freeze-out time tj when the numbers of baryons and antibaryons created are fixed. 
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A rough estimate for 77 may be obtained if we follow the further amplification of the amplitudes Io(titi) of the k = 
modes beyond T\. For large arguments the increase in J (Tm) is mainly due to the exponential exp(rm), so we obtain 



(53) 



For a typical amplification ratio of 5 to 10 during roll-down we then find an approximate freeze-out time of 

U 



This certainly represents a lower limit for the duration of the roll-down, because the increasing <i> 4 contributions to 
the potential will slow down the symmetry-breaking motion. The numerical simulations confirm this simple argument 
and indicate that (mr/) ~ 4—5 provides a reasonably accurate estimate for the freeze-out time (as long as (rom) < 1). 
After the quench, when f§ has assumed its (T = 0)-value f§ — 1, we may neglect the small contributions of the 
explicit symmetry-breaking H to f 2 and to the a-mass m 2 in eqs. JJJ and JHJ and rewrite (|53(1 in the form 

ZL « (54) 
t r m ff 

The typical example for an evolution given in fig|21 shows how during the roll-down the configurations pick up an 
appreciable amount of kinetic energy until the potential starts to deviate from its linear rise and interferes with ((T)). 
Subsequently, ({U}} starts to pick up the unphysical linearly rising lattice contributions l|37(l of the shrinking solitons, 
while the time-averaged ((T)} remains basically constant. As the heavy solitons carry no kinetic energy, ((T)} then 
resides in small transversely propagating fluctuations which eventually are emitted as a and 7r mesons. 



A. Correlation lengths and defect numbers 

In contrast to the integer net-baryon number B — J p d 2 x, the integral (or lattice sum) over the absolute values of 
the local winding density \p\ 

n — J \p\ dxdz (55) 

generally is not integer. The ensemble average of n is closely related to the coherence length R for the field unit- 
vectors in the statistical ensemble of 0(3)-field configurations. If an 0(n)-field is implemented on a (/-dimensional 
cubic lattice with lattice constant a, then the field orientations on the vertices of a sublattice with lattice unit R/a 
can be considered as statistically independent. Then the average ((n)) expected on an N d lattice is 

« n » = v d {aN/R) d (56) 

where v& is the average fraction of the surface of the sphere S d covered by the image of the sublattice unit (this is 
the very definition of a winding density). The number v& can be estimated for different manifolds |17|. For the map 
(compactified-)R 2 — > S 2 defined by the unit-vectors *(x, z) of the 0(3)-field in d = 2 dimensions it is v 2 = 1/4, (i.e. 
l/2 d+1 for each of two triangles which make up each square sublattice unit cell). 

Inserting our result l|33|) (obtained for m = 0) into this estimate for d = 2 dimensions, leads to (( n )) t =t <x T 2 . 
In figO^, this is compared with numbers n and coherence lengths R measured for several initial configurations on an 
N x N lattice for different temperatures (for N — 100 and mass m — 0). Evidently, the finiteness of the lattice causes 
a small systematic deviation from this T 2 -dependence, especially for large values of ((n)), as the coherence length 
approaches the lattice constant. The measured numbers (( n )) follow l|5()|) with satisfactory accuracy for v w 1/5. 
In figOJa the same comparison is shown for non- vanishing mass m = l/£, for l/a = 4. For m 7^ the coherence 
length R can be obtained from i.'>'2\) and compared to the measured values. FiglJb shows that they are reasonably 
well described by R cx T~ 8 . The corresponding measured numbers n follow (|56|) with good accuracy for v m 1/6. 
Of course, small changes of v could be absorbed into a slightly redefined coherence length (note that the correlation 
function jS^I for m — does not decrease exponentially). We shall, however, keep the definition l|3TJ|) . 

The above considerations apply to random configurations which need not contain any fully developed solitons but 
may consist of only small fluctuating local winding densities which cover small fractions of the image sphere. However, 
if the configurations finally have evolved into an ensemble of well-separated solitons or antisolitons embedded in a 
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FIG. 3: Initial (measured) coherence lengths 7? and number of defects n as functions of initial temperature T, measured for 
five random initial configurations for each temperature, on an N — 100 lattice, with m put to zero (left), and with m = 1/t for 
£ = 4 (right). (All quantities in lattice units a). 




FIG. 4: Crosses show the measured coherence lengths R±,R\\ and the average R as obtained from the definition l|34|l . as 
functions of proper time r. For r > 2 they are parametrized by power laws with exponents 0.15, 0.75, and 0.25, respectively. 
The measured defect number n (full line) is compared to the statistical result l|56|l with v = 1/6 (crosses connected by lines). 
(N = 100,T = 0.2,cro = 0.2/ = 4,// = 0.1,A = 1). 



topologically trivial vacuum with only small fluctuations in the local winding density, then the integral Ij55(l counts the 
number of these embedded baryons-plus-antibaryons. We therefore adopt the notion 'number of defects' for (( n )), 
irrespective whether configurations comprise only small local winding densities, or partial or complete solitons. 
For a typical evolution (see e.g. figEJ) the number of defects measured as function of proper time shows a slow 



14 



decrease which follows approximately a power law 

n ~ n (r=To) f^j . (57) 

By the end of the roll-down at freeze-out time 77 this decrease levels off and n settles near the constant which counts 
the number of the finally surviving fully developed solitons-plus-antisolitons (cf. figs^). The decrease in n reflects 
the slow increase in the average coherence length R up till the end of the roll-down. The longitudinal coherence 
length Rn grows very slowly because rapidity gradients are suppressed with 1/r in the Bjorken frame. This leads to 
an effective decoupling of field-vectors in longitudinal direction and subdues the drive for aligning field orientations 
in adjacent rapidity bins. On the other hand, the transverse coherence length R± grows rapidly. For R± 3> R\\ , the 
average radius R obtained from l|34l) is dominated by Ri\. A typical example is shown in fig0]for an evolution which 
starts at To/a — 1. The average R grows with an exponent ofasj 0.25. The statistical argument in (|56() then leads 
to n ~ T~~ 2a , with 2a — 7 ss 0.5. This is slightly steeper than the measured decrease in n. But as the growth in the 
coherence radii sets in only after one or two units of r after the onset of the evolution, the final number of surviving 
defects is reasonably well reproduced by the statistical expression (|56[) (we adopt v = 1/6 from figlHtO- Altogether, 
we typically find exponents 7 ~ 0.4 ± 0.05 for the decrease (|57|) of the number of defects. Then, with l|54|l for the 
typical freeze-out time, we have 

((n)}| T=T/ H<n))| T=Tn (^y' 4 (58) 

for the reduction of the number of defects from its initial value at the onset of the evolution until the end of the 
roll-down. With (to Too-) of the order of 0.5 to 1, we find reduction factors of 1/3 to 1/2, which is not even one order 
of magnitude. So, this is not a dramatic result. The reason is, evidently, that in the expanding Bjorken frame the 
gradient coupling in rapidity direction quickly gets suppressed. 




FIG. 5: The numbers of defects n for two different values H = and H — 1.0 for the strength of the explicit symmetry breaking, 
each for five evolutions on a 100 x 100 lattice (T = 0.3, A = 1, 1 = 5, <7o = 0.2). The arrows point to the freeze-out times t/ f» 22, 
and Tf ~ 11, respectively, where the average lengths |<&| of the chiral field vectors have reached = 1. The initial values of n 
lie within a band from 145 to 185, they all end (for both values of H) in a band from 32 to 52, which corresponds to reduction 
factors of about 1/4. 



It should be noted that all numerically measured exponents are independent of the choice of the lattice constants, 
because scaling x — > ax, (i.e.i — > at), rj — > br\, t — > (a/b)T leaves the EOM i|17|) invariant. The length unit a only 
serves to define the resolution with which the spatial structure of the field configurations is analyzed and all physical 
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results should be independent of this scale. On the other hand, the initial time To denotes the physical point in time 
when the system begins its evolution in terms of hadronic degrees of freedom with a sudden or rapid quench in the 
relevant potential. So, physical results generally will depend on To, as is evident from the reduction factor obtained 
in PI). 

Small explicit symmetry breaking (H ^ 0) accelerates the decrease of n during roll-down, but at the same time it 
reduces the freeze-out time, such that the final number of n remains essentially unaffected by small non-zero values 
of H. FigEl shows a number of evolutions for two different strengths H of explicit symmetry breaking. 

The same is true if additional damping is introduced into the EOM (|17f) by adding a term k8 t <& with damping 
constant n to account for the fact that the field fluctuations are actually emitted from the expanding Bjorken rod, 
carrying away energy. Through this dissipative dynamics the evolutions are slowed down, the roll-down times may 
be retarded by an order of magnitude, but the overall reduction factor in the number of surviving defects remains 
unaffected. Of course, all fluctuations then are damped away during the course of the evolution, and the integrals 
(11311 to l|15l) finally are determined by the remaining ensemble of squeezed solitons alone, while the kinetic energy goes 
to zero. 



B. Meson spectrum 

For times long after the roll-down the average kinetic energy (( T )), and the potential energy parts (( Lj_ + U )) 
(after subtraction of the linearly rising (lattice) contributions from the squeezed solitons l|37f) ) converge towards the 
same constant Ef/2. Their sum Ef represents the average total energy stored in the mesonic field fluctuations after 
the roll-down. Both averages show residual fluctuations around their smooth background with opposite phases, such 
that their sum Ef is smooth. Analysing the spectral density of either (( T )), or (( L± +{/)), (after subtracting the 
background) , tells us about the spectral distribution of pions and a- mesons that will eventually be emitted from the 
expanding Bjorken rod. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 a co 0.9 



FIG. 6: Spectral density n(ui) of the residual fluctuations in the average kinetic energy ({ T )) for times long after the roll-down, 
(N = 80, l/a — 4, H = 0.1). The arrows point to the lowest (double-)frequencies 16H for 7r- and cr-mesons (see text). The 
green line is the exponential exp(— 12au>). 

We consider the Fourier-transforms 

c(w) + is(u) = J T \{ T(r) - T(r) ))e^dr, (59) 

where the integral covers times long after the roll-down, e.g. r a /ro ~ 100, t^/tq ~ 1000, and T(t) subtracts the 
smooth background. The absolute value, e(uj) — V c 2 + s 2 , represents a spectral energy density, from which we may 



16 



extract the spectral particle number density 

n(w) = e M = 53 „W f( w - 2u;«) % - 2u#>) + • • , (60) 

The r/-sum with z, j = 0, 1, 2, ...N/2 covers all frequencies on the lattice for pions and cr-mesons with masses m T and 
m a given in (jHJ 

" ^(^) 2 + (^) a+ ™^ ■ <«> 

Generically, T(r) contains contributions ~ (cos(w|jV)) 2 from the pionic fluctuations, and ~ (cos(lu^ t) + c) 2 from the 
cr-fluctuations around some nonvanishing average c. Therefore, the spectral functions e(uj) and n(uj) will, in addition 
to the double frequencies 2iJf^ , also contain contributions for the cr-mesons at the single frequencies Jf^ . 

Figure HO shows the spectral density n(u) as obtained from the residual fluctuations in the average kinetic energy. 
The long vertical arrows point to the first four 2wjj pionic frequencies (|61|l for ij = 00, 10, 20, 30, with rn^ — H/£ 2 , 
(H = 0.1, l/a = 4, /o = 1). It may be seen that the overwhelming part of the strength resides in the lowest and 
first excited pionic modes. The strength decreases rapidly with excitation energy, approximately like exp(— 12aw). 
The same is true for the strength of the cr-modes. (The short arrows in figEl point to the first three modes with 
ij = 00,10,20, with single frequencies iJf^ and double frequencies 2w|j ). However, the number density ^2 n ij ' f° r 

the sigmas, which we may extract from the strength located at the double frequencies 2a; jj', is only about 5% of the 
pionic strength residing in the first three pionic modes. For an order-of-magnitude estimate of the pionic multiplicities 
we therefore ignore the c-contributions. 



C. Meson and baryon multiplicities 

To obtain a simple estimate for the energy Ef finally available for meson production we consider the time of the 
onset of the roll-down t\ in (|51|l which marks the transition from the gradient-dominated to the potential-dominated 
phase. At this time, for sufficiently small Uq, the total energy is dominated by the linearly rising term (r /tq)CqVo in 
the potential l|44|l . With the onset of the roll-down the average potential ({{/)) starts to deviate from this linear rise 
and bends down to interfere with (( T )) and (( L )), (cf. fig. |2J). In the numerical simulations the large-time limit of 
(( U )) and (( L )) is masked by the (lattice-artificial) rise of the soliton contributions. But the asymptotic (( T )) 
is free of these artifacts and (apart from residual fluctuations) approaches a constant value, which is well represented 
by the linearly rising \{t / tq)CqVq taken at r = r%. Approximating n by 5/(v / 2w (T ) as given in (|51|l . we then have 
(with f 2 = 1) 

E f = /^C Vo« ff^%(abN% (62) 
to 8y2 

where we have again neglected the small contribution of the explicit symmetry-breaking H to the tr-mass. Within this 
level of accuracy we can also ignore that about 30% of the pions carry the energy (instead of m w ), and obtain 
the pion multiplicity n T from dividing (|62|l by m^, 

"'-fl{abN 2 ). (63) 



8^2^ 



This number may be compared with the baryon-plus-antibaryon multiplicity given in l)58[l. We use for (( n ))i T=To 
the statistical result l|56[). with initial spatial coherence length i?o- Then we have 

<<n>> = ^YfeY. (64) 



Ro J V 4^2 

The last factor relies on the estimates ()51f) and (|53fl for the times T\ and 77, which are valid as long as (Tom) < 1, 
(otherwise they have to be obtained more accurately from l|50|l and l|52|ll. The evolutions described above have been 
performed for initial configurations selected with netbaryon number B — 0. So, the average number rip of antibaryons 
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created during the phase transition is (( n ))/2. With typical values v ~ 1/4. -) ~ 0.4 we find for the multiplicity 
ratio of antibaryons to pions 



a m„ (t to ct ) ■ 

lp /n^0.U- — JI . (65) 



With an overall energy scale / 2 of the order of the pion mass , and Rq of the order of m 1 = \[2m a 1 , this ratio is 

np/n^ w 0.07(^m cr )(r m (T ) 7 . (66) 

The ratio a/6 of the spatial and rapidity lattice constants which appears in this result has a physical meaning: 
according to (|31|) it is equal to the ratio of the (transverse) spatial coherence length Rq and the (longitudinal) rapidity 
coherence distance Ri\q in the initial configuration. Naturally, this ratio is of the order of tq. So, for initial times To 
typically of the order of the inverse er-mass we find antibaryon-to-pion multiplicity ratios of the order of 0.05 to 0.1. 



VI. GENERALIZATION TO 3-D 0(4) 

For the generalization to the 3+1 dimensional 0(4)-model we keep the parametrization as given in eqs.Q, ||3J|, 
and Jfjjl. In this case /J is an overall constant of dimension [mass 2 ], so the physical fields /„■<!? are of mass-dimension 
one. The winding density is no longer given by but we keep as defined by the second equality in eq.l|^)l. 
Conventionally, the strength of the £( 4 )-term in © is given in terms of the Skyrme parameter e as 

2 A/ 2 1 

(67) 



(8^) 2 4e 2 / 2 ■ 

In this case the typical spatial radius of a stable skyrmion in its rest frame is mainly determined by the balance 
between and £' 4 ', so it is of the order of (e^) -1 . 

For the map (compactified-)R 3 — > S 3 defined by the unit vectors of the 0(4)-field in 3 spatial dimensions the 
statistical result <|56|) for the average number of defects found on a (aN) 3 lattice for initial configurations with 
coherence length Ro generalizes as 

««»|r=r = ^ T (aiV/ J Ro) 3 . (68) 

(We again use a = br® for the lattice constants). The factor 5 counts the number of 3-simplices (tetraeders) which 
make up a cubic sub-lattice cell of size Rq, the factor (l/2 d+1 ) with d = 3 is the (absolute value of the) average 
surface area covered by the image of one 3-simplex on the image sphere S 3 . So the factor 5/16 counts the average 
'number of defects' associated with a cubic lattice cell with lattice constant given by the initial coherence length Rq. 
A certain arbitrariness in the definition of the coherence length may translate into modifications of this factor 5/16; 
(e.g. for a random lattice of 3-simplices the factor 5 is replaced by (247r 2 /35 ~ 6.8) I n anv case we do not 

expect order-of-magnitude changes in this factor as compared to the d — 2 case, where we had 2/2 2+1 . 

However, through the cubic power the result is now very sensitive to the actual value of Rq in the initial ensemble. 
Different concepts about the physical nature of the initial configurations will imply quite different ways to arrive at the 
appropriate initial coherence lengths Rq. For an initial ensemble which is characterized by a temperature T we could 
proceed as in and relate Rn to the temperature, or to the mass m 2 (T) = X\f 2 (T)\/£ 2 of the field fluctuations; 
but it has also been suggested [22j to tie i?o to the parton density which makes it independent of the temperature 
concept. So, for the moment it seems appropriate to keep the initial coherence length Rq as a parameter. 

Adding a second transverse dimension does not change the result (|34|l for the average of the transverse and longi- 
tudinal coherence lengths. The growth in the resulting R again is dominated by the slow increase of the longitudial 
coherence length i?y which is unaffected by additional transverse dimensions. The estimates (|51|) and <|53|) for the 
times T\ and r/ of the onset and end of the roll-down also remain unaffected, as they only rely on the amplitudes A(t) 
of the transverse waves (|48f) . irrespective of the number of spatial dimensions. (In this case, 

£(4) 

now contributes to 

L± with a term containing four transverse gradients, acting on the direction of the 0(4) field. The roll-down, however, 
takes place in areas which are topologically trivial, i.e. with small angular gradients, so we do not expect a strong 
effect on the roll-down times. Within the approximations which led to Q58[l. we then find for the average number of 
baryons+antibaryons present after the roll-down 
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with a ~ 0.2 to 0.25. We denote the transverse area (aN) 2 of the Bjorken rod by A, and replace the ratio (a/b) of the 
lattice constants again by the initial time r . Then we obtain for the rapidity density of antiprotons (n p — j(( n ))) 



drip 5 tqA ( rom CT \ 
~drj = 64 ~Rj \AV2 J 



3a 

(70) 



(At this point we count all baryons as nucleons, assuming that excitational fluctuations and rotations contribute to 
the surrounding pionic fluctuations). Although the strength of the Skyrme term does not appear explicitly in l|70|l . 
the presence of the £( 4 )-term is essential for the formation of the solitons, because during the evolution it transforms 
the sum over the absolute values of average winding densities into the average number of fully developed solitons and 
antisolitons. So it is hidden in the growth law characterized by a. 

For meson production we adopt the considerations which led to the estimate <|62ll . Counting again all mesons as 
zero momentum pions we have 

«x = ^= ^CoVo = f^^(AbN). (71) 
Wtt m n t 8V2 m n 

The rapidity density of negatively charged pions n n - = \ri- K then is 

drj '" r 24^m 7r v ' 

In a heavy-ion collision the transverse area A of the Bjorken rod will correspond to the spatial overlap of the colliding 
relativistic nuclear slabs. As we have assumed spatially homogeneous initial conditions we have to consider slabs with 
constant nucleon (area-)density. In order to account for the number A of nucleons contained in one slab, its radius 
must be taken as r^A 1 / 2 , with ro ~ 1.2 fm. Then, as function of centrality, (dn^- /drf) is directly proportional to 
the number of participants N p , which is one of the basic experimental results in relativistic heavy- ion collisions. For 
central collisions of A-nucleon slabs we have N p = 2A, so we find for the 7r"-rapidity density per N p /2 participants 

1 dn w - 5n m a t r 2 



N p /2 drj 24V2m n 



(roUy. (73) 



This is an interesting result because all parameters have been absorbed into physical quantities. There are, however, 
several caveats: We have used for this result the form of the potential © after the quench (only this enters into 
the calculations). This means that differences between the average potential energy before the quench l|39[) and 
immediately after the quench (|44l) are left out. However, this difference is of the order a 2 , which has been neglected 
in (|73|l anyway. But generally, the result l|73l) should be considered as a lower limit. It should further be noted that 
the result l|73|l depends linearly on the time t\ for the onset of the roll-down. The definition of t\ in (|51|) is not 
very stringent and may be subject to changes by ±20%. The estimate (cf. eq.lJSJ)) we used for the time T\ required 
(toTO ct ) < y/2 , i.e. with m a ^3 — 5 fm -1 , the initial time To should not exceed 0.3 — 0.5 fm. Another unsatisfactory 
feature of the homogeneous Bjorken rod is that the inhomogeneity in the nucleon (area-)density of real relativistic 
nuclear slabs with transverse radii roA 1 / 3 has to be represented through radii roA 1 / 2 for homogeneous slabs. 

The experimental value |l6j | for the 7r _ -rapidity density per N p /2 lies between 1.2 and 1.5 for N p increasing up to 
350. With rp/jr = 0.57, m a /m v 5 — 8, our result l(75|) leads to 

' ' / "' 0.75- 1.2. (74) 



N p /2 drj 

In the light of the reservations discussed above, this is quite satisfactory. 

In contrast to the parameter-free pion multiplicities, the result for baryon-antibaryon creation depends on two 
parameters: the time To when the initial hadronic field ensemble is established and begins its expansion, and the 
initial coherence length R within that ensemble. From (|70|l and l|72|) we have (with a ~ 0.2 to 0.25) 

-g- « 0.15 y ° } . 75 

The experimental value for the ratio of integrated p to n~ multiplicities lies between 0.065 and 0.085 ^6] for varying 
numbers of participants. With m^/f 2 = 3.0 fm, and a typical c-mass of m a m 3 fm -1 , the experimentally observed 
multiplicity ratios are reproduced if Rq and To (both in [fm]) satisfy 

Ro w (3t ) q+1 / 3 . (76) 
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For initial times in the range 0.2 < tq < 0.5 the dependence on a is very weak and the coherence length varies in 
the range 0.7 < Rq < 1.2 (all in [fm]). These values are certainly within the limits of conventional assumptions. 
Interpreted in terms of a thermodynamic equilibrium ensemble, Rq ~ lfm ~ T _1 implies the standard estimate 
T ~ 200 McV for the chiral phase transition. With our choice a/b = tq = Rq/ R^q for the ratio of the spatial and 
rapidity lattice constants, the initial time tq ps 1/3 fm resulting from l|76l) for Rq = 1 fm then means that in the initial 
ensemble the initial rapidity coherence distance Rm extends over three units of rapidity. 

VII. CONCLUSION 

We have presented numerical simulations of the dynamical evolution which chiral field configurations undergo in 
a rapidly expanding spatial volume. Starting at an initial time To from a random hadronic field ensemble with 
restored chiral symmetry, we follow its ordering process and roll-down into the global potential minimum with spon- 
taneously broken chiral symmetry. In accordance with standard concepts of heavy-ion physics we have considered 
one-dimensional longitudinal expansion of an essentially baryon-free region of high energy density, as it may be realized 
in the aftermath of an ultra-relativistic collision of heavy ions for central rapidities. 

Performed on a space-rapidity lattice in proper time of comoving frames, such simulations are very powerful in- 
struments which allow to investigate a multitude of interesting features related to the chiral phase transition. We 
have concentrated here on the topological aspects which are directly related to baryon-antibaryon multiplicity as a 
sensitive signal for the phase transition. Mesonic abundancies could be analysed as well, both for tt and a mesons (or 
any other elementary fluctuations included in the chiral field). Not only their spectra can be obtained, but from the 
instantaneous configurations the spectral power of their momentum distribution could be extracted at every point in 
time. 

The method is not restricted to thermally equilibrated initial ensembles with global or local temperature; in- 
homogenities and anisotropy in the correlation lengths could be implemented naturally. Surface effects could be 
investigated by suitable boundary conditions. This may be interesting with respect to the A-dependence of spectra 
and multiplicities. Here we have applied only standard periodic conditions. The onc-dimensional expansion could be 
replaced by anisotropic or spherically symmetric expansion, which may be of specific interest in cosmological appli- 
cations. We have used the sudden quench approximation, which could be replaced by any desired time-dependence 
of the chiral potential with arbitrary quench times. We have selected ensembles with conserved net-baryon number 
B = or very small B. Any other choice would be possible, and it appears as a peculiarly attractive feature to study 
evolutions in ensembles with high net-baryon density, either fixed or in the form of grandcanonical ensembles. The 
method is well suited to analyse distribution, growth and realigning of domains with disoriented chiral condensate 
as has been shown previously in purely dissipative dynamics [l4|. The generalization to SU(3)-fields appears most 
interesting, to learn about strangeness production in terms of baryonic and kaonic abundance ratios. 

Evidently, the method opens up a wide field of applications. Unfortunately, however, we know very little about the 
nature and characteristics of the initial ensemble which enters crucially into all physical results. So, in our present 
analysis of antibaryon and pion multiplicities, the experimental data do not allow to draw definite conclusions about 
the validity of the topological approach, because the results depend on two initial coherence lengths, the spatial Rq 
and rapidity R\\q, (which for an isotropic initial ensemble are related by Rq = R\\qTq). We can only conclude that 
conventional assumptions about these quantities lead to results which are compatible with experimentally detected 
multiplicities. So, luckily, the mechanism is not ruled out. On the other hand, an assumption like tq = Rq ( which 
would imply that the correlations have grown with the speed of light from a pointlike origin ) is ruled out: it would 
overestimate the abundance ratio in (|75() by a factor of 5. 
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